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ABSTRACT 



O 

^ ■ We compute the canonical partition function for quantum black holes in the approach 

bi) of Loop Quantum Gravity (LQG). We argue that any quantum theory of gravity in which 

the horizon area is built of non-interacting constituents cannot yield qualitative corrections 
to the Bekenstein- Hawking (B-H) area law, but corrections to the area law can arise as 
a consequence additional constraints inducing interactions between the constituents. In 
LQG this is implemented by requiring spherical horizons. The canonical approach for LQG 
CN [ seemingly favours a logarithmic correction to the B-H law with a coefficient of — |. Our 

CN I initial calculation of the partition function uses certain approximations that, we show, do 

^ \ not qualitatively affect the expression for the black hole entropy. We later discuss the 

quantitative corrections to these results when the simplifying approximations are relaxed 
and the full LQG spectrum is dealt with. We show how these corrections can be recovered 
^ ■ to all orders in perturbation theory. However, the convergence properties of the perturbative 

■ series remains unknown for now. 
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I. INTRODUCTION 



Classical and semi-classical properties of black holes have been studied quite vigorously 
in past several decades As a result, we have discovered that some nice thermodynamic 
properties can be associated with them. For example, black holes are now understood to 
be endowed with a temperature and to possess an entropy characterized by their mass 
and by any other charges that may be associated with them. These properties, although 
derived on the semi-classical level, are expected to be robust enough that explaining them 
in terms of an ensemble of microstates has come to be considered of prime importance for 
any quantum theory of gravity. Many approaches to quantum gravity, like String theory [Q , 
Loop Quantum Gravity (LQG) and the theory of Causal Sets f^, have therefore taken 
up this problem and each has succeeded in some measure in providing some deeper insight 
into the problem. However, the Bekenstein-Hawking area law seems generic enough that 
all the approaches, although differing considerably in design, are able to reproduce it. As 
a result, there remains no general consensus on a final theory of quantum gravity and it 
becomes essential to attempt not only to obtain the original semi-classical results, but also 
to discover, if possible, signature corrections of quantum origin to the semi-classical wisdom. 

In this paper, we take up one of the popular approaches to quantum gravity, namely LQG. 
Black holes have been studied quite extensively in this approach f^, mostly in context of 
isolated horizon scheme 0. In LQG, the area spectrum of an isolated horizon is readily 
available and black holes with a given horizon area have been studied in the microcanonical 
ensemble [[r|-pi]|]. Logarithmic corrections to the Bekenstein-Hawking entropy have been 
shown to exist on the semi-classical level |Tl[| and as a signature of LQG 0, |10[- They 
have been proposed in other approaches as well [jl2| but there is some confusion about the 
precise value of the coefficient of this term. Approaches that employ conformal symmetry 
techniques generally seem to prefer the value —3/2 whereas the pure quantum geometry 
approach seems to favor the value of —1/2. 

We analyze LQG black holes in the canonical area ensemble with a finite number of 
punctures, A^. The use of the area ensemble was first advocated in [|l^ and has since 
been employed by others (see, for eg. IH, W^ and references therein). A justification for 
its use has recently been proposed in [|1^. The use of A^ as a statistical variable in the 
area canonical ensemble was first advocated in [jl^ and its importance for black holes was 
emphasized in There are similarities between our approach and that of |TD| and our 
canonical calculation agrees with the results therein, although the latter work was performed 
in the microcanonical ensemble. The microcanonical and canonical entropies do not have 
to coincide, of course, except in the thermodynamic limit. Now it is well known that 
the microcanonical entropy in LQG suffers from the fact that the black hole entropy is 
not a differentiable function of its arguments but rather a ladder (or staircase) function 
1^ , p!9| pO[] . As pointed out in this makes it difficult to interpret basic thermodynamic 



variables, such as the temperature, which are defined in terms of derivatives of the entropy. 
The canonical partition function is assumed to be, by contrast, a smooth function of its 
arguments. We obtain the canonical partition function without making any assumptions 
concerning our variables or referring to any semi-classical feature, so the following analysis 
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is independent of any d priori thermodynamic input. We expect that our approach will 
provide a useful framework for examining area fluctuations about the microcanonical value 
and for a generalization to the grand canonical ensemble, where we should be able to precisely 
test some of the assumptions of [|1^, |1^. Here, we find general agreement with |1T0[ . 



In section II we briefly review the LQG approach to black holes. We begin our compu- 
tation of the canonical partition function in section III by making a certain approximation, 
which we will refer to as the "shell" approximation and define later, to the LQG area 
spectrum. Ignoring the projection constraint, we easily recover the Bekenstein-Hawking 
area-entropy relation together with an estimate of the Barbero-Immirizi parameter and we 
show that any correction to the Bekenstein-Hawking entropy is inherently absent. In sec- 
tion IV we impose the sphericity (or projection) constraint on the black hole horizon and 
re-evaluate the partition function. The effect of the projection constraint is to produce a 
logarithmic correction to the standard semi-classical result. Section V is devoted to a dis- 
cussion on relaxing the "shell" approximation and incorporating the full LQG spectrum in 
a perturbative fashion. We develop a systematic approach to compute the corrections to all 
orders. We conclude by discussing some issues raised as well as the general outlook related 
to our approach in section VI. 



II. BLACK HOLES IN LQG 

As a non-perturbative, background independent approach to canonical quantum gravity. 
Loop Quantum Gravity (LQG) has met many challenges effectively, although a full dynam- 
ical picture of quantum gravity has remained elusive so far. Starting with holonomies, built 
from su(2) valued connections, and fluxes, built from densitized triads, as basic variables 
of the theory, one obtains a well-defined Hilbert space made up of cylindrical functions 
acting on spin-networks over which the holonomies and densitized triads are defined. A 
spin-network is a graph with edges labeled with su(2) representations and nodes character- 
ized by su(2) intertwiners. Spin-networks are eigenstates of the Area and Volume operators 
and this leads to a reasonably good description of some of the geometrical observables of 
the theory, at least on the kinematical level. For example, an edge with spin-representation 
j carries an area of eigenvalue 

A, = Stt^II^jU + I), (1) 

where j G and 7 is an unknown parameter of quantization known as Barbero- 

Immirizi parameter. Alternatively, an area element of area Aj can be thought as being 
punctured by an edge of the spin- network carrying representation j. In general a surface 
that is punctured by many edges of different representations will have the area spectrum 

A = 87r7/2^Vj,(j, + l), (2) 

i 

where the sum is over all intersections of the edges with the surface. One considers the black 



hole horizon as an isolated horizon [21] which is threaded by many edges of a spin-network. 
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FIG. 1: Horizon punctured by the edges of a spin-network 



Each edge carries some representation and hence deposits some area to the horizon as shown 
in figure |1]. 

In the microcanonical ensemble, one simply tries to count the number of different con- 
figurations, i.e., the different ways of threading a horizon with a fixed area A by spin- 
representations. This then gives the black hole entropy. However, counting the total number 
of microstates in this way has not been trivial owing to the fact that the area spectrum in 
(0) is non-distributive in nature. This issue has been discussed and treated in different ways 
by various authors 0-0. In this way, one obtains the asymptotic formula 



where we have put the Barbero-Immirzi parameter explicitly in the area as it appears in the 
expression for the area, and utilizes the flexibility of flxing 7 to recover the Bekenstein- 
Hawking result. 

In the following two sections we will examine LQG black holes in the canonical ensemble. 
We begin by making the "shell" approximation at flrst, but later, in section V, we generalize 
to the full LQG spectrum. This leads to corrections to the entropy derived in sections HI and 
IV but these corrections amount to a renormalization of the Barbero-Immirzi parameter and 
do not change the functional dependence of the canonical partition function on the extensive 
variables. 

III. SETUP FOR THE CANONICAL CALCULATION 

Let us begin by flrst ignoring the projection constraint and, since all punctures are con- 
sidered to be distinguishable, let us assume that there are rij punctures carrying spin j, 
then 



AA(7) 



(3) 





(4) 
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with ttj capturing the information about the area spectrum. For a configuration containing 
distinguishable punctures, rij of which carry spin j, we will have 



distinct configurations. The degeneracy factor (2j + 1)"^ signifies 2j + 1 values for rrij = 
{—j, —j + 1, .., j} for a spin- representation j. Therefore, we may write the partition function 

as 

(5) 



where /3 is conjugate to the area and is not the black hole temperature. By the Binomial 
theorem, this is 

N 

(6) 



Now, since j G {|, 1, |, . . .} it is convenient to make the change 2j = I so that 



\l=l,2,... J 



N 



0- 



(7) 



where a = Air'jlpP. We will now study the general properties of this kind of partition 
function without restricting ourselves to to this or any particular area spectrum. 

A. Formal Solution 

If, as we have above, 

Zia,N) = z'^ia) (8) 
then, because A = it follows that 

A din z{a) dcf 

= q. (9) 



An-fllN da 
Inverting this relation, one finds a solution of the form 

a = 47r7/2/3 = a(g). (10) 
The entropy of the system is obtained from 

^ = In Z{a, N) + /3A{a, N) = N[\n z{a) + /3a(a)] (11) 
where, (3 = f3{q) as above and A = Na i.e., a = in'ylpq, as can be seen from (|^). Thus, 

S{q) = N[lnz{a{q)} + qa{q)]. (12) 
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Now let us maximize the entropy with respect to the total number of constituents (punctures 
for the LQG model), A^. Then 

where the prime refers to a derivative with respect to q. Now from (|^), 

dq A q 



ON 47r7/2iV2 



(14) 



so, 



S 



o + a + ga 



and, since 



we have 



or 



da 

d\iiz A 



da ATT-fl^N 
S 



: 0, (15) 

-q, (16) 



^ Q^iQ) = (17) 



In z{a{q)} = ^ z{a{q)} = 1. (18) 

We must solve this equation for N = N{A) and reinsert into S to determine the entropy as 
a function of area. In this way we will get 

S= N{A)qa{q)\^^^^^^^^,, (19) 

but the relation z{a{q)} = 1 implies that q = qo, a. constant. This gives 

^(^) = A \ > (20) 

and 

A 



The entropy is always proportional to the area so long as the canonical partition function 
obeys (^) and regardless of the spectrum. It follows that any theory of quantum gravity in 
which non-interacting constituents make up the horizon (like the punctures in LQG at this 
level) will result in the Bekenstein-Hawking area-entropy law. To obtain the Bekenstein- 
Hawking relation in LQG framework, the value of the Immirizi parameter should be fixed 
at 

7 = (22) 

TT 

Within this framework, any correction (in particular logarithmic) must arise from additional 
constraints (such as the projection constraint in LQG) on the horizon. Therefore, we do not 
expect any logarithmic correction to Bekenstein-Hawking relation in shell collapse scenarios 
||22|| . Now we estimate the numerical value for 7 for an approximated LQG spectrum, which 
looks very much like the spectrum in dust shell collapse. 
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B. The 'Shell" approximation to the LQG spectrum 



We define the "shell" approximation by y/P + 21 = ^y (I + ly — 1 / + We have called 
this the "shell" approximation because an identical solution is obtained in the canonical 
quantization of the LTB dust models [^. In those models the collapsing dust ball is viewed 



as made up of dust shells, and the apparent horizon for each shell is shown to have a similar 
area spectrum. In making this approximation, we will be making errors towards low spin 
punctures, but for large spin punctures, the errors will not be that significant. Then, 

= ^ (; + l)e-<'«) = e-««) = (23) 

1=1,2,... 1=0 ^ ) 

SO that the condition z{a{q)} = 1 translates to 

2 - e"" = (e" - 1)^ (24) 
There are two possible solutions for e^, viz., 

and the equation for q 

dlnz ^ 2 1 

determines go in terms of these solutions. There is only one positive solution for g and it 
corresponds to e°" = 2.247; we end up with 

a(go) = 0.810, go = 3.318, 7 = 0.258, (27) 

using (P2|). It is important to note that even if an exact calculation could be done then the 



entropy would still be proportional to A while the constants go and 7 would change. Now, 
for the estimate of the Barbero-Immirizi parameter from microcanonical counting, Ghosh 
and Mitra, in their study of LQG black holes [Q, suggest that the correct value for 7 is 
0.274. On the other hand. Ling and Zhang produce a very similar value, 7 ~ 0.247, from 
their study of = 1 super-symmetric LQG black holes. Importantly, this value of the 
parameter 7 is closely related to the agreement between LQG spectrum and quasi-normal 
(ringing-) mode frequency of black holes 0. One can in principle, get a better estimate for 
7 by dropping the "shell" approximation. We will illustrate how this can be done in section 
V. 

For the present we recall that the projection constraint is still missing and in the next 



section we will introduce it. This constraint has been identified as the possible source 
of a logarithmic correction from analyses in the microcanonical ensemble. First we write 
the canonical partition function with this constraint, but still staying within the "shell" 
approximation. 
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IV. THE PROJECTION CONSTRAINT 



We now consider a spherically symmetric black hole. Such a system in LQG will be 
described by a spin-network which is an eigenstate of the projection operator fi ■ J with zero 
eigenvalue. In other words, the nij values carried by the spins puncturing the horizon must 
add up to zero. This is the projection constraint. It can be written as, 



X]^im,"^i = 0, (28) 



where njmj gives the number of punctures carrying spin j and projection = — j, —j + 
1, . . . , j — 1, j along the axis. We can see that the punctures are now "interacting" in 
the sense they have to satisfy a constraint for getting the desired quantum state and their 
interaction will manifest itself in terms of some correction to 
First we write, 

S^njmj.mj = p, (29) 

for some integer then we will put the constraint on p, forcing it to vanish. Suppose there 
are = Ylj mj ''^jm.j punctures, then the partition function becomes 



Z(/3, N)= J2 n I ^P'' e-'^''^^-^ (30) 



where SpQ is the Kronecker delta function. We want a suitable representation for this 
function. For integer p a convenient representation isQ 

^^^,0 = 7^ rdke^'^P. (31) 

and using it we find 

Z(/3, AT) = — V ( r rffce'*'^-. "^™^'"^\) e-'""^^^^-^ "-"^'^^ (32) 

Interchanging the sum over rijm and the integral over k, 



re 



To see this, note that the right hand side is 

1 e'''P 



2tt ip 



1 [g2.ip_l] 



2'Kip 



Since p € Z, if p ^ then the right hand side is vanishing. In the hmit as p — > 0, this becomes 



hm = 1, 

p->o 27rp 



so the representation is true. 
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N 

2lT 



i- r dkij2 e(2^^"^^ -«-^^"^ ) , (33) 



again after using the Binomial theorem. Note that to recover the result with no projection 
constraint, we must drop the integration over k and set k = 0. Now let 2j = I as before, 
2mj = ri & {—I, — / + 2, . . . , / — 2, /} and consider the sum in the integrand first. 



°° Jk{l+2) _ -ikl 



E ( E ) = E ' (34) 

1=1 \ n / 1=1 
We easily check that the limit as A; — of the sum is 



lim ' ... = / + 1, (35) 



and therefore ( |5^ ) reduces to (|^. 

We want to get the partition function 



N 



N 

ik{l+3) _ -ik{l+l) 



^ (3g) 



into a manageable form. In the second step, we made use of shell approximation to simplify 
the calculation. After a little manipulation, we end up with the expression 



1 r / 2 cos A; - e-'^ ^ ^ 



Z(I3, N) = — dk{ , (37) 

where we have made use of the evenness of the integrand to rewrite the limits of the in- 
tegration. Note that this is not of the form (^ and that this is because of the projection 
constraint, which introduced an interaction between punctures as discussed above. In the 
limit as /c — 0, the integrand approaches 

7^;r^^ (38) 



which is precisely ( p3D as required. 

Now we try to approximate the integral from the information we have about the integrand. 
First, we see that it is a unimodal symmetric distribution between its limits. So we would 
like to approximate it by a "bell-curve" . (In order to do so, we can make the transformation 

A; = 2tan-^(x/2), (39) 
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to get the limits form — oo to oo.) We then make an assumption that it is well- approximated 
by a Gaussian in some regime of A^. If true, then the variance (cr^) of the "bell-curve" is 
obtained by the second derivative of the integrand at the peak, since in the case of a Gaussian 



In our case. 



1 / 2 cosk(x) — e 



2tt ye^'^ - 26"^ cos k{x) + 1 
1 2 cos k{x) — e 



TV 



dk 



dx 



27r(l + xV4) Ve^'^ - 26" cos k{x) + 1 

resulting in the variance 



N 



x2 



cr 



/(x)|o 2(6^^- 1)2(26'^ -1) 



f"{x)\o -1 + e'^(4 + e'^(-5 + e'^(2 + AN))) ' 
which in large N limit tends to 



(41) 



(42) 



2 e'^-12 2-e-'^ , 
lim ^ 43 

Now, being a symmetric distribution its skewness is already zero. We analyze its kurtosis 
next. For a Gaussian with a zero mean, the kurtosis obtained from fourth standardized 
moment, can be given by its distribution function as 

. /i4 f{x)y^'Kx)\o . .... 

[/"(x)loF ^''^ 

which gives the "excess kurtosis" as 

/32 - 3 = 0. (45) 
For our case, also due to the fact that its mean is zero, we obtain, 

/M _ /(x)|o/W(x)|o _ 6[(1 - 2e")^(e" - 1)^ + 8e^"(-l + e"(2 + e'^je" - 1)))N + Se^^N^] 
^ ~ [f"{x)\o? ~ [-1 + e-(4 + e-(-5 + e-(2 + 4A^)))]2 

(46) 

and we see that in limit of large A^, the excess kurtosis of this distribution also vanishes, 

hm ^-3^0. (47) 

For large iV the integrand tends to a unimodal symmetric distribution with zero skewness 
and vanishing excess kurtosis. Thus, our initial assumption is justified and in the large iV 
limit this function is indeed well approximated by a Gaussian. In this limit the partition 
function is given by the area of a gaussian, which can be readily evaluated. Therefore, we 
have the partition function in the thermodynamic limit 

Z(iV,a) = /(a;)|oV^a. (48) 
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or, explicitly, 



and thus, 



Z{N,a) 



y - 1)2(2 -e-^) / 2-e 



(e- - 1)^ 



N 



In Z In z{a) — In + const., 
(assuming that is large) where z{a) is given in (0). We find, as before, 

A d\n.z def 



A-K-tllN da 



and 



Now 



S = \nZ + f3A = N[ln z{a) + crq] - - In + const.. 



^•^1 1 1 n 
-— : = Inz — ^ Inz = (J. 

ON 2N 



(49) 



(50) 



(51) 



(52) 



(53) 



for large A^ and retaining terms up to the order of 1/ yiV. This implies, once again, that 
z{a) = 1, and it can be solved for go and o"(go) as before. Then 



S ^ qoa{qo)N{A) - - In N{A) + const. 



(54) 



with A^ = A/A'K'yllqQ, i.e., 



47r7/2 2 V4vr7/2go 
or with the chosen value of the Immirizi parameter 



const., 



(55) 



(56) 



Thus we obtain a logarithmic correction, whose origin lies clearly in the imposition of the 
projection constraint. The logarithmic correction comes with a factor — i, as already pointed 



out in 1^, 10 1, and therefore the canonical calculations are shown agree with the microcanon- 
ical results. The microcanonical ensemble and the canonical ensemble agree even at the 
subdominant correction, in contrast with the results of [|I^, in which the coefficient is de- 
termined to be +|. The reason for this discrepancy is that the authors in [jl5| work in the 
grand canonical ensemble (with a vanishing chemical potential). Fluctuations in the number 
of punctures can be shown to contribute precisely a logarithmic term with coefficient +1 
to the entropy. The sum of this contribution and the logarithmic term from the projection 
constraint in the canonical ensemble above leads to their result. The projection constraint, 
absent in the shell-picture, plays a pivotal role in bringing about the LQG signature. 
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V. TOWARDS A FULL z(cr) 



We will now try to see how the full LQG spectrum is likely to improve upon the calculational 
details. In the previous sections, we made the approximation 



v/^(7TT) = + 1)2-1 + 



(57) 



in evaluating the partition function. As already argued, this type of approximation will 
induce a significant error only towards low spin punctures. For a more precise calculation, 
we actually want 

^(^) = E + l)e--v^(^^, (58) 



1=1,2,. 



which can be re-written as 



[a) = Y: {l)e-^\ 



(59) 



1=2,. 



The most direct way to evaluate the sum above is to employ the Mellin-Barnes representation 
of the exponential function. 



/ ds T(s) a-' 

2m 



(60) 



where r G M"*", Re(Q;) > and the integral is taken over a line parallel to the imaginary axis. 
The integration path can be closed in the left half plane. 
Using this representation we find 



z[a)= y — I ds Vis) a-n-'+^ 

l_2 Jt-ioo 



v 



1=2 



ds T(s) a~n~'+^ 



1 + 



fe=0 



2{k+i) s{s + 2)...{s + 2k) 
2 X 4 ■ . . . X 2(A: + 1) 



) 2'n-i 
1=2 
oo ^ 

=E— 

1=2 ^r-ioo 
1 pT-\-ioo 

'^li J T~lOO 

^ 1 /"T+ioo 

+ E ^ / r(s) a-^iCis + 2A; + 1) - 1} 



fc=0 



s{s + 2)...{s + 2k) 



2 X 4 ■ . . . X 2(A; + 1) 



(61) 



In the above equations ({s) is the Riemann zeta-function and T{s) is the Gamma function. 
Now we may systematically expand, order by order in k, to estimate the possible corrections. 



A. Shells again 

It is easy to verify that the first sum yields the "shell" approximation, for it can be 
written in the form 

z{a) ^ — / ds T{s) a-^{as - 1) - 1}, (62) 
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which has simple poles ats = {2,0,— 1,...}, so has the value 



n=0 



using residue theorem. Now, using the relation 

a-n) = - 



Bn+1 



n + 1 



between Riemann zeta-function and the Bernoulli numbers, we get 



0"^ ^ — ' n\ [ n + 2 

n=0 



and again, starting with the expression 



n=0 



cr — ' n! 

n=2 



we find 



and therefore 



E 

n=0 



(n + 2)n! Scr ^1 - e^'^ 



z{a] ^ -— ( ^ I - e = ^ + . . . 



da \1- e-" 
This is precisely what we have been using. 



(e- - 1)^ 



B. First correction 



To the above we must add the infinite series in k. For /c = we have two terms 
/ ri.s .srr.s^ rr-^r.s + 1^ - ^ / 

47rz 



— I ds sT{s) a-X{s + l) - — ds sT{s) a 

^'^'^ J T—ioo 'iTTl J r—ioo 



The first term has a simple pole at s = with residue +1 and also simple poles at s 
Z~ with residues 



(n- 1)! 



so, it IS 



n=l 



(n- 1)! 



(-1)"B„<T' 



n=l 



cr 



2(1-6-: 
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using (|6^) and (^). The second of these is also straightforward since the poles are only 
those of sr(s), which lie at {—1, —2, . . . , — n, . . .} with residues {—l)"'~^/{n — 1)! so 



fT+lOO 

Ani 



1 



n=l ^ ' 



^ (-l)V" a 
2^ ^ = ^e' 



n=0 



Combining the two, this correction is 



a 
2 



1-6-" + e 



-2a 



and so 



2 - e"'" a 
ie- -1)2 + 2 



1 - e-'^ 

1 - e"'" + e~2<^ 



1 - e-'^ 



(72) 

(73) 
(74) 



The condition z{(y) = 1 becomes a transcendental equation and cannot be solved in closed 
form but one can numerically estimate the left hand side and we find only one positive 
solution for go, 

a(go) ~ 1-187, go ~ 0.510, 7 ^ 0.378. (75) 

This gives o"(go) and hence 7 values that are somewhat higher than in the previous order. 
Obviously, the series oscillates about its asymptotic value so an important question to ask 
is: how fast does the series converge? We will deal with convergence in a subsequent work. 



C. k>0 

Each of the next correction terms would be of the form 

"^^'^ ^ s{s + 2) ...{s + 2k) 



1 pT-\~tQO 

— dsT{s) a-'{C{s + 2k + l)-l} 



(76) 



2 X 4 ■ . . . X 2(A; + 1) 
with A; G N. For example, at order /c = 1, the correction term would be 

— - / ds s{s + 2)T{s) a-Xis + 3) - — - / ds sis + 2)r(s) a-^ (77) 

which can be evaluated as before. The first integral has simple poles at s = —n e 
{ — 1, —2, . . .}. with residues 

(_l)n^(n-2)C(3-n)a" 

j , « ^ 2, (78) 

n! 

and residue —a^ at n = 2. The second integral has simple poles at s = —n G { — 1, —2, . . .} 
with residues 

n\ ' 

and so on and so forth. All the higher order corrections can be computed and subsequently 
added following similar calculations. 

It is hoped that the larger the order, the better will be the estimate, but this would have 
to be rigorously shown from a convergence analysis. This issue is not taken in this paper, 
but will be discussed separately elsewhere. 
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VI. CONCLUSION AND DISCUSSION 



In this work, we have analyzed quantum black holes in the framework of the canonical 
ensemble. We first obtained the Bekenstein Hawking relation and demonstrated that this 
relation is a generic feature of any quantum gravity theory built around non-interacting 
constituents of horizon area. Such theories will not show any significant departure from the 
semi-classical analysis, apart from any numerical estimates of parameters involved in the 
quantization scheme. We also studied black holes in the Loop Quantization scheme under 
some simplifying approximations (the "shell" spectrum) whose effects are expected to be 
smeared in case of large black holes. For those black holes we determined the Barbero- 
Immirizi parameter to be quite close to one found in literature, relating it to the ringing 
mode of black holes. 

For the particular case of LQG, it was shown unambiguously that the theory predicts a 
logarithmic correction to the Bekenstein-hawking entropy and that this correction is a direct 
consequence of spherical horizons. Our computation leads to a logarithmic correction with 
a factor of —1/2 and not —3/2 as has been suggested in the literature [|1^. The reason for 
this is that we are working with the reduced U{1) Chern-Simons theory for which there is 
one projection constraint. The factor of —3/2 is obtained in models which employ the full 
SU{2) symmetry in which there are three constraints [Q. The factor of —1/2 agrees with 
the results of g [ig. 

The partition function for the full LQG area spectrum was obtained as an infinite series 
of terms, making use of the Mellin-Barnes representation of the exponential function, and 
corrections at various orders were discussed. Convergence of the series remains an open 
question for now and will be tackled in a future publication. Interestingly, this approximation 
scheme was avoided in with the use of Pell equation. The exact calculation will help 
understand the convergence of the series and the full and exact subleading character of the 
entropy. We are presently examining this possibility and will report on our results elsewhere. 
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